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Регуляризация дифференцирования 
с использованием репагулярного 
вейвлет-преобразования 


В статье рассматривается проблема регуляризации дифференцирования изображений с помощью 
репагулярного вейвлет-преобразования, для решения которой выполняется представление этого 
преобразования с помощью интегрирования дробного порядка. 


Введение 


Одной из базовых процедур в системе компьютерного распознавания зритель- 
ных образов является сегментация изображений на однородные по какому-либо 
признаку или набору признаков области (сегменты). Целью этой процедуры является 
уменьшение объема обрабатываемой при распознавании информации за счет выделения 
границ однородных областей. Составной частью процедуры сегментации изображений 
является подчеркивающее преобразование, задача которого — концентрация энергии 
сигнала вблизи границ сегментов. Именно подчеркивающее преобразование определяет 
главным образом помехоустойчивость процедуры сегментации и погрешность выделе- 
ния границ однородных областей. Кроме того, при распознавании иерархических 
объектов, т.е. имеющих структуру «объект-подобъект», а также при анализе сцен под- 
черкивающее преобразование позволяет выделять объекты требуемых геометрических 
размеров. Для таких задач требования к помехоустойчивости процедуры сегментации и 
погрешности выделения границ объектов зависят от уровня иерархии последних. 

В составе преобразования, подчеркивающего перепады значений признака сег- 
ментации, часто используется дифференцирование. Задача дифференцирования изобра- 
жений является некорректной [1], [2]. Чтобы избежать проблем, связанных с некоррект- 
ностью задачи дифференцирования, последнее выполняют в пространстве обобщенных 
функций. Для этого производят свертку изображения, подлежащего дифференцирова- 
нию, с гладкой основной функцией, которая дифференцируется перед выполнением 
свертки. Такая последовательность действий реализуется несколькими концептуально 
различными способами. 

1. Строится гауссовское масштабно-пространственное представление изображе- 
ния, Т.е. однопараметрическое семейство сглаженных изображений, линии уровня кото- 
рых являются выпуклыми. Деталям изображений большого масштаба соответствуют 
детали изображений малого масштаба, причем это соответствие устанавливается с по- 
мощью гауссиана. Гладкость гауссиана позволяет определить в смысле обобщенных 
функций производные изображения, для которого построено масштабно-простран- 
ственное представление. 

2. Используется регуляризация по Тихонову, т.е. ищется дифференцируемая функ- 
ция, близкая в смысле некоторого критерия к изображению признака сегментации. 
Минимизация соответствующего критерию функционала производится с помощью 
регуляризации, учитывающей порядок производной. 
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3. Применяется линейная фильтрация, реализующая дифференцирование, для 
подчеркивания перепадов значений признака сегментации. Обычно предполагается, что 
задана модель перепада значений признака сегментации, а также критерий качества 
подчеркивания перепада [3]. 

В работе [4] выведена взаимосвязь между дифференцированием с помощью ре- 
гуляризации по Тихонову и гауссовским масштабно-пространственным представле- 
нием, а также показано, что фильтры Канни-Дерича могут применяться для решения 
задачи регуляризации дифференцирования. 

В случае иерархической структуры объекта распознавания для получения мас- 
штабно-пространственного представления в [5] применялось репагулярное вейвлет- 
преобразование (ВП). Это преобразование определено в [6] как свертка строк или 


|“ зепх, [| <=,; 


столбцов изображения с функцией и(х,а) = где =, — фиксиро- 


0, х>=,, 
ванное число, зависящее от ае (0,1) — параметра преобразования. Параметр а репа- 
гулярного ВИ определяет соотношение между помехоустойчивостью и погрешностью 
выделения перепадов значений признака сегментации. Применение репагулярного ВП 
для подчеркивания перепадов значений признака сегментации позволило достичь 
более высокой помехоустойчивости по сравнению с применением дифференцирова- 
ния [6]. Поэтому целесообразным является представление результата репагулярного ВП 
как регуляризованного решения задачи дифференцирования. 

Математический аппарат, используемый в данной работе, — операции дифферен- 
цирования и интегрирования дробного порядка, которые представляют собой более 
общий подход к подчеркиванию перепадов значений признака сегментации, включаю- 
щий и обычное дифференцирование. Эти операции использовались наряду с обычным 
дифференцированием в ряде работ [7], [8] для подчеркивания перепадов интенсив- 
ности и других задач обработки изображений. 

Целью данной работы является доказательство того, что применение репагу- 
лярного ВП осуществляет регуляризацию дифференцирования. Для этого выполняется 
представление репагулярного ВП с помощью интегрирования дробного порядка. 


Представление репагулярного вейвлет-преобразования 
с помощью интегрирования дробного порядка 


Определение интегралов и производных дробного порядка зависит от свойств 
функции, стоящей под знаком производной или интеграла, а также от области определе- 
ния этой функции. В качестве области определения функции, стоящей под знаком 
производной или интеграла, рассматривают конечный отрезок, полуось или ось. 
Предполагают, что функция, стоящая под знаком дробной производной или интегра- 
ла, является интегрируемой по Лебегу или обобщенной. 

В данной работе дифференцирование или интегрирование дробного порядка при- 
меняется к функции 2(х), представляющей значения признака сегментации (напри- 
мер, интенсивности) строки или столбца изображения. Предположим, что эта функция 
задана на конечном отрезке [а,Б] и может быть нулями дополнена на всю вещественную 


ось так, чтобы интегралы, участвующие в определении дробных производных, сходи- 
лись на бесконечности и в точках отрезка [а,5] . Например, выберем &(х)е Г,(К), где 


1< р<!/а, а - порядок дробной производной или интеграла. 
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Дробные производные на всей вещественной оси вводятся в [9] следующим обра- 
зом. Выражение 


ь 1 Я Е ®(6а 
ре, 
Га-а) 4х +, (х-й 
называется левосторонней дробной т. порядка а, а выражение 
1 (ща 
(р“=)(®) = о 
Га-а) 4х > (х-0 
— правосторонней дробной производной порядка а, 0<а<1. 
Для > 1 полагают 


ре 
око 


< х<о 


по 


[повода м [@] +1. 


Ги-а) 4х" у 
Дробные интегралы в случае и. оси обозначаются как 
(118) =—— | ВОР. хо. 
Г = 3 (х-й“ 
г аа 


(в 8) о ее ©<х<о, 


и называются интегралами дробного порядка & > 0, соответственно, левосторонним 
и правосторонним. 

Эти интегралы определены, например, на функциях 8(х)е [,(-0°,°) при 0<@<1 
и 1<р<1!/а. 

Представим репагулярное ВП в терминах дробного интегрирования по всей ве- 
щественной оси. 

те Га ра Г Г) 
(“= сье) = Г 8(1) -] 8) +| 8(1) 


х-й| зеп(х-Ю 5 х-Й зеп(х- И) х— |" ет(х- 0) 


— ыы 
т (= Им Ре - геи $) ©) Га) (Г $) ©. 


В частотной области репагулярное ВП в терминах дробного интегрирования 
записывается с помощью формул преобразования Фурье дробного интеграла: 


Е(1-“ 8) (6) = (©) (16) * , 0<1-а<1, (1) 
Е(Г“ 8) (6) = 86) 6)", 0<1-а<1, (2) 
А @ 
где (+10)* = |" е2“ ‚ В — вещественное число, которое в данном случае равно |1-@, 


Е - оператор преобразования Фурье. 
С учетом формул (1), (2) репагулярное ВП в частотной области имеет вид: 


те а -а © © 
Е(аткдона))=Р(е{и в) Пед(рк)о) Геи Г, 
Обозначим результат репагулярного ВП функции 2(х) через Ё(х): 
Е(х) = (“р х,о). 
Тогда формулу (3) в частотной области можно Пе. в виде: 


Яо) = пево —т =) 


Со“ 6“ 


. (3) 
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Преобразуем (е) т. о., чтобы в сумме осталось только одно слагаемое: 
—(1-а) т 


а) а -(-а) а 
Ее) = |@ е ? -|@ е ? - ®(@)Г(&)= 


ПН а а 


К 
5-Я зале 27 оугде 2 Е = (4) 


—(1-а)л 


= (10) ^®. 6(0).е 2 "т о о.) | 
Учитывая формулы (1), (2), получаем: 
©) «деи 38 - о Тур наи И о Иа 
= ГЕ (198). (©), (5) 
где # (@)=2щ? ва оо (-одл зап @). 


В пространственной области 


№ (х) = 25° 29 я 


9(х)-2 


81((1-а)л) 
- : 


Регуляризованное решение задачи дифференцирования 
изображения с применением репагулярного вейвлет- 
преобразования 


В [1 показано, что задача дифференцирования функции &(х), известной прибли- 
женно, является существенно некорректной. Для решения существенно некорректных 
задач применяются методы регуляризации [1]. 

Покажем, что регуляризованное решение задачи дифференцирования функции 
2(>х) ‚ известной приближенно, может быть получено в результате применения репа- 
гулярного ВП. В процессе этой операции производится дифференцирование регуляри- 
зованного сигнала изображения. 

Использование регуляризации в случае численного дифференцирования заклю- 
чается в интерполяции и аппроксимации данных некоторой функцией с последующим 
вычислением производных этой функции [2]. Предположим, что интерполяционная 
или аппроксимирующая функция и(х) для значений строки или столбца изображе- 


ния такова, что и’(х) = Е(х). Тогда (Р'и)(х) = /(х). Следовательно, 


Ка) = о). (6) 
Используя соотношения (4) — (6), сформулируем следующую теорему. 
Теорема. Если функционал 


Е = | (- (о) 8) + (1) Ков (во, 


где параметр ^.>0 достигает экстремума при @)=&(0), то в пространственной 
области ч,(х) может быть представлена как линейная свертка функции 2(х) е Г (К) 
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с импульсной характеристикой фильтра Н(х), которой в частотной области соответ- 
ствует передаточная функция 


Несут, 


—(1-а) 7 
по 1- 
где (о) =2е ? т” т ша о, 


Доказательство. Из представления (5) следует, что 


Ио) = ГЕ (6) (6) = Г(едсюу  ® во) (6). (7) 
Подставив в (7) соотношение (6), получаем 
[еее Кое) =0. (8) 
Г(@а) 


где Ко) =! (6). 
Уравнение (8) представляет собой условие равенства нулю первой вариации 
функционала 


Е) = Сбобо оС’ Воде оао, ©) 


Т.К. 


Е 


ае 8(6)=в (6) — |8 т Г(@) (озер 


где # (6) - функция, на которой функционал (9) достигает экстремума. Последняя 


#(о)=& (6) › 


может быть получена путем умножения преобразования Фурье &(@) исходной функ- 
ции 2(х) на #(6) = Ге) В (9). 


Обозначим Я = . Заметим, что А > 0, в случае @ > 0. Тогда (9) принимает вид: 


2Г (а) 


Е) = | (= (<) 8) + 1-6)“ Ко) (ао, (10) 


а #(@)= Роу оо) . Теорема доказана. 


Функционал (10) является результатом формирования целевой функции задачи 
оптимизации в области преобразования Фурье. Покажем, что этот функционал прини- 
мает вещественные значения. По аналогии с доказательством предложения 1 работы [4] 
предположим, что результат &(@) преобразования Фурье исходной функции ©(х) 
является вещественным, а функция ©(х) симметрична относительно начала координат. 
Так, если е(х) представляет значения признака сегментации строки или столбца изо- 
бражения, ее можно продолжить симметрично относительно начала координат. При 
таком предположении функция /(х) как результат репагулярного ВП 2(х) антисиммет- 


рична относительно начала координат. Следовательно, преобразование Фурье (ео) функ- 
ции }(х) является чисто мнимым. Тогда согласно (6) функция (@) принимает вещест- 


венные значения. Выражение (—1)'” <) также является чисто мнимым потому, что 
(—160)"" (—10)—“ 
= - = — $ 
р 1-а)л . (@-а)л 
(6) ег о о, 25ш М ап 


1-а 
й 


(-ю)"“ Ко) = 
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Тогда выражение (—)“ Ио) вещественно, т.к. 


о" 


25т аи ——_—^_ 5210 
2 
Следовательно, функционал (10) принимает вещественные значения и его опти- 
мизация проводится в пространстве вещественных чисел. 


В процессе доказательства теоремы было получено выражение для функционала 
в частотной области, экстремум которого достигается на функции & (4). Результатом 


(-16)“ Во) = (-6)С-ю)"* К) =-@ 


обратного преобразования Фурье от &(4) является функция и,(х) в пространственной 


области. Получим выражение для функционала в пространственной области, экстремум 
которого достигается на функции и, (х). 


Согласно равенству Парсеваля [10] для функций 4(х), ъ(х) и результатов пре- 
образования Фурье этих функций имеем: 


Гасдуок=5 | пожемо, (и) 
ее ты 


где черта над функцией означает комплексно-сопряженную величину. Равенство Парсе- 
валя выполняется при условии, что 4(х), у(х), &{@), @) интегрируемы с квадратом 


на оси К. 
Применим равенство Парсеваля к (10). Сначала предположим, что &@) = &(@), 


а ©) =). Так как ©) — вещественная функция, то 6) = (0) ‚ а 4(х) = ®(х) и 
у(х) =и(х). Следовательно, 


[ &е)Щедае = 2 | а(х)у(х)ах =2л | 2(х)и(х)ах. (12) 


Далее применим (11) ко 2-му слагаемому в (10). Полагая &(@) = 4(—10)”“ ® (&),а 
©) = &), имеем а(х) = д(6:=” и). (х), е)=Йо). 

Получим представление функции <) в пространственной области. Для этого 
преобразуем Ио) следующим образом. 


.(1-а) п ЗИ 
Ио) = ® (@) -е ы т"! |“ и вп о, . = 2 кв“ ии еп О] . (13) 
1 


(1-а)л 


Тогда й(х) = ты 108 2 


С учетом (13) функция Же) в пространственной области имеет вид: 


Ро | 2 +в ша В] в 
2. 2 х 2 
Используя (14) для 2-го слагаемого в (10), получаем 
| ее 0)“ в (о) Ко)ав = м 
т 1 1-@/2 ›@ ал . 
Е о” и) 62969 -4е ие ах. 
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Подставляя (в и (15) в (10), имеем 


Еи(х)]=-2л | и(х)е (ах + 2л4 [2 рее и). 62969 ь се - я Р я ах = 
Хх 


—© 


=-2л | исдасдах+дяя | (р-*Зи } (х)0(х)ах- о м -8)% | [(5: т и). - т 


—с 


+14, | [Ра ) (х)ах. 


С учетом известного соотношения [10] | у(х)о(х)ах = (0) ‚ где у(х) — некоторая 


функция, Е[и(х)] принимает вид: 


Е[и(х)] = -2=| и(х)®(х)ах +474. [-ы ) (0) -8л4с1ю их | (5: и. и). в Ча + 
+4лА, [ (ри ) (х)ах. 


— 


Введем параметр /, = 21Й м, где и-— чисто мнимое число, такое, что Ди и < 0. То- 
гда А 20, причем А Е А. Видоизмененный с учетом значения А, функционал Ё[и(х)] 
обозначим Ё[и(х)]. Он имеет вид: 


Е [и(х)]=-2л | и(х)е (хх +2) л | (реги) (хх+4йл (ри) (0)— 
_ Е. (16) 
(-а)л [ (2: 2/2 и) ©) 1. 


—© 


—874 с 


Интерпретация подчеркивающего преобразования 
как результата регуляризации задачи 
дифференцирования изображения 


Методы регуляризации в задаче сегментации изображения используются с целью 
избежать некорректности решения, вызванной наличием шума. Это происходит, когда 
область значений подчеркивающего преобразования изображения не замкнута или 
когда задача является существенно некорректной. Последнее означает, что изменения 
правой части операторного уравнения задачи обработки изображения, связанные с на- 
личием шума, могут выходить за пределы образа множества исходных изображений 
при отображении его упомянутым оператором. Так как изображение неизбежно содер- 
жит шум, методы регуляризации являются хорошим инструментом для решения задачи 
сегментации изображений. 

Основная идея метода регуляризации Тихонова заключается в том, что вводится 
последовательность непрерывных приближений не непрерывного оператора. Регуляри- 
зация обобщенного решения осуществляется через однопараметрическое семейство 
непрерывных операторов К, , ^ >0, таких, что при ^,->0 К, &(х) стремится к обоб- 


щенному решению. Применение К, к незашумленному изображению дает приближен- 
ное решение, которое улучшается при 7, —> 0 с точки зрения качества приближения. 


Однако если изображение содержит шум, имеем К, в, (х) = К, 8(х)+ К, п,(х), где 
о 
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2,.(х)=а(х)+п,(х), п,(х) - шум. Слагаемое №, п,(х), отвечающее за распростра- 


нение ошибки приближения, расходится при 4-0. Поэтому вводится функционал 
Ф(>.(х))-+ ЛЧ (2, (х)) и метод регуляризации путем выбора параметра Л, выполняет 
согласование между качеством приближения, определяемым слагаемым Ф(з, (х)), и 
распространением ошибки, выражаемым слагаемым (2, (х)). 

Доказанное в этой работе предложение утверждает, что проинтегрированный 
результат репагулярного ВП исходного изображения доставляет минимум функциона- 
лу (16), представляющему целевую функцию оптимизационной задачи на безусловный 
экстремум. Представим эту задачу на безусловный экстремум как результат применения 
метода неопределенных множителей Лагранжа к вариационной задаче на условный 
экстремум. Тогда последняя реализует вариационный принцип отбора возможных ре- 
шений для метода регуляризации, применяемого к некорректной задаче дифференци- 
рования изображений. Сформулируем вариационную задачу на условный экстремум, 
для которой функционал (16) представляет собой функцию Лагранжа. Для этого за- 


[ее] 
метим, что минимизация скалярного произведения —2л | и(х)®(х)ах по и(х) может 
проводиться путем минимизации квадрата евклидового расстояния между этими 
функциями. Известно, что репагулярные вейвлет-функции характеризуются разной 
регулярностью, охватывая диапазон от производной дельта-функции б’(х) до функции 


Хевисайда 0(х) [11]. Путем свертки с дб’(х) осуществляется обработка изображения 
в дифференциальных методах контурной сегментации, тогда как @(х) является мо- 


делью перепада интенсивности для корреляционно-экстремальных методов. 
В данной работе предполагалось, что исходная функция ©(х) представляет строку 


или столбец изображения, а /(х) — результат репагулярного ВП функции ©(х). Тогда 
перепадам значений признака сегментации функции ©(х) соответствуют пики /(х), 
ширина которых зависит от значения & , параметра репагулярного ВП. 

Функция и(х) определяет проинтегрированный результат репагулярного ВП. 
Поэтому перепадам значений признака сегментации функции 2(х) соответствуют сгла- 
женные перепады значений признака сегментации и(х). Чем шире пики /(х), тем более 


[ее] 
сглаженной является функция и(х) и выше значение —2л | и(х)=(х)ах . В этом случае 
—© 


результат репагулярного ВП приближается к результату согласованной фильтрации. 


Следовательно, слагаемое —2л | и(х)=(х)ах при минимизации функционала (16) опре- 
“> 
деляет качество приближения функции ©(х) функцией и(х). Последнее влияет на 
такую характеристику подчеркивающего преобразования, как погрешность определения 
координат точек перепадов значений признака сегментации изображения. 
Известно, что функция Лагранжа строится для оптимизационной задачи с огра- 
ничениями в виде равенств. В процессе определения вида этих ограничений, прежде 


1-а)л инь о 1 
всего, слагаемому 8лй, сам | (5! “?и) (х)—@х поставим в соответствие огра- 
х 


—о 


ничение [ (ри ) м =0: 
№ 


— 
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Это условие означает, что квадрат левосторонней дробной производной от про- 
интегрированного результата подчеркивающего преобразования должен затухать на 
бесконечности, убывать в окрестности нуля, а также должен быть симметричным 
относительно начала координат. Тогда, с точки зрения обработки изображения, пик- 
сели, расположенные по обе стороны контура, равноценны. 

Условие убывания квадрата левосторонней дробной производной от проинте- 
грированного результата подчеркивающего преобразования в окрестности нуля позво- 
ляет, приравняв 3-е слагаемое функционала (16) к нулю, сформулировать ограниче- 


ние (27“?и)* (0)=0. 


ее) 
г и 
Значение слагаемого 24л | (5! Ри (х)4х возрастает при уменьшении 
—® 


ширины пиков результата репагулярного ВП функции ©(х). Тогда это слагаемое 
определяет распространение ошибки приближения функции 2(х) функцией и(х) и 
влияет на помехоустойчивость подчеркивающего преобразования. 


Слагаемое 2) л | [9* 2); (х)4х в выражении (16) представляет собой стаби- 


лизирующий функционал при вариационном подходе к регуляризации задачи дробного 


2 
дифференцирования изображений. Это обусловлено тем, что [Ри (х)>0 и >0 


тт 


как неопределенный множитель Лагранжа. Вариационный принцип отбора возможных 
решений задачи регуляризации осуществляется путем минимизации именно стабилизи- 
рующего функционала, который интерпретируется как условие гладкости проинте- 
грированного результата подчеркивающего преобразования вдоль контуров объектов 
на изображении. Физическим основанием для этого условия является то, что в силу 
свойств процесса обработки незашумленному изображению соответствует ограничен- 
ный спектр и, следовательно, оно имеет ограниченную производную [2]. 


Выводы 


Получено представление репагулярного ВП с помощью интегрирования дробного 
порядка. Доказано, что применение репагулярного ВП осуществляет регуляризацию 
дифференцирования изображения. Следовательно, рассматриваемое преобразование 
может использоваться для подчеркивания границ однородных областей при построении 
методов сегментации изображений в задачах распознавания иерархических объектов, а 
также при анализе сцен. Направлением дальнейших исследований предполагается 
построение масштабно-пространственного представления изображений с применением 
репагулярного ВП. 
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М.В. Полякова 

Регуляризащя диференциюовання з використанням репагулярного вейвлет-перетворення 

У статт! розглядаеться проблема регуляризаци диференциювання зображень за допомогою репагулярного 
вейвлет-перетворення, для вирипення яко! виконуеться представлення цього перетворення за допомогою 
1нтегрування дробового порядку. 
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